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Abstract: Generalized Word Length Pattern (GWLP) is an important and widely- 
used tool for comparing fractional factorial designs. We consider qualitative factors, 
and we code their levels using the roots of the unity. We write the GWLP of a 
fraction P using the polynomial indicator function, whose coefficients encode many 
properties of the fraction. We show that the coefficient of a simple or interaction 
term can be written using the counts of its levels. This apparently simple remark 
leads to major consequence, including a convolution formula for the counts. We 
also show that the mean aberration of a term over the permutation of its levels 
provides a connection with the variance of the level counts. Moreover, using mean 
aberrations for symmetric s'" designs with s prime, we derive a new formula for 
computing the GWLP of P. It is computationally easy, does not use complex 
numbers and also provides a clear way to interpret the GWLP. As case studies, we 
consider non-isomorphic orthogonal arrays that have the same GWLP. The different 
distributions of the mean aberrations suggest that they could be used as a further 
tool to discriminate between fractions. 
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1.1 Introduction 


In design of experiments. Generalized Word-Length Pattern (GWLP) is an important tool for 
comparing fractional factorial designs. For a regular fraction of a full-factorial design P with 


m factors, the Word Length Pattern (WLP) of P has been introduced by Suen et al. (19971 as 
the sequence A{P) = {Ai{P), A 2 {P), ■ ■ ■, Am{P)), where Aj is the number of defining words 
with length j. Such a measure of the degree of aliasing can be easily interpreted in the regular 
case. The GWLP has been generalized for non-regular asymmetrical designs by |Xu and Wu| 


(2001 but it has a less evident meaning than in the regular case. 


The aberration and the GWLP through the polynomial indicator function of the fraction 


P have been introduced in Li et al. (20031 and Cheng and Ye (20041 for two- and three-level 


cases respectively. In those papers the aberration of a simple or interaction term is defined as 
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the square of the module of the corresponding coefficient of the indicator function, and the j-th 
element Aj (7^) of the GWLP is the sum of the aberrations of the terms of order j, j = 1,... .m. 

As demonstrated in Xu and Wu (20011, the GWLP does not depend on the choice of a 
particular orthonormal basis of the functions defined over T>, while the aberration does. |Pistone| 


and Rogantin (20081 use the complex coding of the factor levels to express the basis of the 


functions, and in particular of the indicator function. With this coding the coefficients of the 
indicator function are related in a simple manner to many interesting properties of the fraction 
and allows us to define aberration and GWLP in a clear way. The complex coding is particularly 
useful in the case of qualitative factors, as assumed in this work. To simplify the computation, 
avoiding the use of complex numbers, Fontana and Pistone (20131 represent the coefficients using 
the counts of the levels appearing in each simple or interaction term. As general references for 


GWLP and its properties, the reader can refer to Mukerjee and Wu (20061 and Ghen and Cheng 

pol^ . 

The practical use of the GWLP to discriminate among different designs is well known. 
Given two designs J-i and T 2 , the Generalized Minimum Aberration (GMA) criterion consists 
in the sequential minimization of the GWLP. is better than T 2 if there exists j such that 
A\{Ti) = A\{T 2 ), ■. ■ ,Aj{Ti) = Aj{T 2 ) and Aj+i{Ti) < Ajj^\{T 2 ). Despite the fact that 
the GMA criterion is widely applied, the statistical meaning of the elements of the GWLP is 


somewhat unclear. In the original work of Xu and Wu (20011, the GWLP of a symmetrical 
design is written as the MacWilliams transform of the distance distribution. This result has 
been generalized in |Qin and Ai| ( [2007[ ) to the case of multilevel designs. Under a different point 
of view, Gromping and Xu (20141 writes the first non-zero element of the GWLP as the sum 
of the B? coefficients of suitably defined linear models. The connection between GWLP and 
Discrete Discrepancy has been investigated in |Qin and Fang (20041. 

In this work we use the expression of the GWLP via the aberrations of the interaction 
terms of a given order. In turn, the aberrations are computed using only the levels counts of 
the corresponding terms. We fully exploit such new expressions in two directions. First, we 
establish a convolution formula for the counts of the terms, in symmetrical s’” designs with s a 
prime number. Second, we introduce the mean aberration of a simple or interaction term, over 
the permutations of its levels. The mean aberration has a very simple expression, and it is easy 
to compute and to explain. Indeed, we show that the mean aberration is proportional to the 
variance of the level counts. 

Moreover, we prove that for symmetrical s’” designs, s prime, the j-th element Aj(iF) of 
GWLP is the sum of the mean aberrations of the terms of order j, j = 1,... ,m and therefore 
the mean aberrations produce an alternative decomposition of the GWLP. In our knowledge, 
the proposed formula is the simplest over all alternative expressions in literature. Nevertheless, 
in general, this property does not hold. 

The paper is organized as follows, fn Section [l.2| a short review of the algebraic theory of 
factorial designs is given. In Section [l.3| the convolution formula that expresses the relationships 
among the level counts is obtained, while in Section |1.4| the mean aberration is defined, and 




































1.2. ALGEBRAIC CHARACTERIZATION OF FRACTIONAL DESIGNS 


its properties and its connection with the GWLP are studied. Section |1.5| is devoted to the 
comparison of fractions that have the same GWLP but different distributions of the mean 
aberrations. Finally in Section [l.6| we briefly describe some directions for future work. 

1.2 Algebraic characterization of fractional designs 

In this section, for ease in reference, we present some relevant results of the algebraic theory of 
fractional designs. The interested reader can find further information, including the proofs of 
the propositions, in [Fontana et al.] | |2000[ ) and jPistone and Rogantin] | |2008| ). 

Let us consider an experiment which includes m factors. 

Let us code the Sj levels of the j-th factor by the s^-th roots of the unity ^ = 

exp ^ k^, k = 0,Sj — 1, j = 1,..., m. We denote such a factor by Qsj , = 

{wo, • . . , 

As a = P mod s implies ui^ = , it is useful to introduce the residue class ring and 

the notation [fe]s for the residue of k mod s. For integer a, we obtain ( 0 ;*,)“ = uj[ak]s- We also 
have oJhOJk = We drop the sub-s notation when there is no ambiguity. 

We denote by D the full factorial design with complex coding: 

T> = Di X ■ ■ - Vj ■ ■ ■ X Vm with Vj = Qsj ■ 

the cardinality of the full factorial design is = YYjLi ■Sf • 

We denote by L the exponent set of the complex coded design {0 ,..., Sj — 1}, j = 1,... ,m: 

L = X • • • X Zs„ . 

Notice that L is both the exponent set of the complex coded design and the integer coded 
design. The elements of L are denoted by a, /3,...: 

L = {a = (qi, ... ,a,n) : Oj = 0,..., Sj - 1, j = 1,... ,m} ; 

[a — P] is the m-tuple ^[qi — ,..., [aj — Pj]^. ,..., [om — . The computation of the 

j-th element is in the ring Z^^. 

In order to use polynomials to represent all the functions defined over T>, including counting 
functions, we define 

• Xj, the j-ih component function, which maps a point = (C, • • •, Cm) of to its j’-th 
component, 

Xj : T> B (Cl,..., Cm) ' ^ Cj T’j . 

The function Xj is a simple term or, by abuse of terminology, a factor. 

• X°‘ = Xf‘^ ■ . .. ■ Xfp", a £ L = Zsj X • ■ • X Zs^ i.e., the monomial function 

A“ (Ci,...,Cm) ^ cr-.-.-c™ • 


The function A“ is an interaction term. 
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The following proposition provides the level set of a term for any choice Si, ..., Sm of 
the levels of the factors. 

Proposition 1. On the full factorial design T>, the simple term takes values in Qtj, where 
tj = Sj/ gcd(aj, Sj). The interaction term takes values in where ta = lcm(ti,..., tm) 
and ti (i = 1,... ,m) are determined as before. 

Proof. It follows from the properties of the arithmetic in □ 

We observe that for a symmetric s"* design with s prime number, t = s for all a £ L. 

The set of monomials {X°‘ : a G L} is an orthonormal basis of all the complex functions 
defined over P. 

Definition 1. A fraction is a multiset {iFt,R) whose underlying set of elements J-* is con¬ 
tained in P and R is the multiplicity function i? : —>■ N that for each element in R, gives the 
number of times it belongs to the multiset R. We call R counting function. 

The underlying set of elements J-t, referred to as the support of Rt, is the subset of P 
that contains all the elements of P that appear in T at least once. We denote the number of 
elements of a fraction T by n, with n — ^(0- 

We use the basis : a € T} to represent the counting function of a fraction. 

Definition 2. The counting function i? of a fraction is a complex polynomial dehned over 
P so that for each ( £P, i?(C) equals the number of appearances of f in the fraction. A 0 — 1 
valued counting function is called indicator function of a single replicate fraction R. We denote 
by Ca the coefficients of the representation of i? on using the monomial basis {X“, a £ L}: 

RiO = C&R,cc.£C. 

aeL 

Proposition 2. Let T he a fraction with counting function R. 

1. The coefficients Ca of R are given by: 

^ I] uJh (1.1) 

^ C6.F ^ h^O 

where ta is determined by X°‘ according to Proposition^ and Ua.h is the number of the 
occurrences of Wh in {X“((C) : ^ £ R}. 

In particular co = n/ffP. 

2. If R is a single replicate fraction, the coefficients Ca satisfy the following relationships: 

^ ( Cg C\^a — g] (1'2) 

geL 


for each a £ L. 



1.3. CONVOLUTION FORMULA WITH COUNTS 


Proof. 


1. We have 




C6.F 


#V 


#V 


TT'Ct.ta—h 


where the hrst equality is proved in Pistone and Rogantin (20081, the second equality is 


proved in Fontana and Pistone (20131, and the last one derives from properties of the 
roots of the unity. 


2. See Pistone and Rogantin (20081. 


□ 


The coefficients Ca encode many interesting properties of the fraction T as orthogonality 


among factors and interactions, regularity, and aberration, see Pistone and Rogantin (20081. 


1.3 Convolution formula with counts 


It is well known that the simple and interaction terms are dependent each other. The proposition 
below gives such relationships in term of the level counts of all terms on a fraction T. 

Proposition 3. Let T be a single replicate fraction of a symmetric designs P, with s prime 
number. The counts Ua, Ua = (na.o,..., Ua.s-i), a,re related according to 


s-l / s-l \ 

^ ' I ^ ^ ^ ^ —/3],[fc —*] 1 0 

fe=0 \ PeL i=0 j 

for each a £ L. For convenience we define no,fc = 0 for k = 1,..., s — 1 


Proof. From Equation (1.2 I we have 


C[s-c,] = ^ C[s-p] C[^_c] . 
peL 

Substituting the coefficients C[s_q,] expressed in terms of the level counts (ua.o, 
in Equation ( |1.1| , we have: 

s —1 ^ s —1 s —1 

fe=0 " PSL i=0 j=0 

Taking k = i+ j, the thesis follows. 


Corollary 1. We denote by ra.k the coefficients of uik in Equation (1.3l; 

S —1 

'f'(y.,k — ^ ^ ^ ^ ‘^13,i '^[a — l3],[k — i] ■ 

i=Q 

If s is prime, for each a, Proposition^ gives s — 1 relationships among counts: 

^a,0 ~ = • • • = rct,s — l • 


(1.3) 


i,s —1)5 


□ 
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Proof. If s is prime, a polynomial X]fc=o '''kOJk is zero if and only if ro = ri = • ■ • = fs-i- □ 

We observe that previous relationships are not independent both because of fraction prop¬ 
erties and complex numbers properties; for instance; 

C2 n(o,,,.,o),o = n and n(o,,.,,o),fe = 0, for fc = 1,..., t - 1; 

C3 ^£12 ,... ,a^), A; [^i — ai ], [s 2 — ( 12 ],..., [st,i — ), [s^ — A;] fot (ui, 0.2., • • • 5 Om) ^ {0; • • ■ 5 Si 1} X 

{0,..., S 2 — 1} X • • • X {0,..., Sm — 1} and k £ {0,..., t — 1}. 

C4 possible relationships coming from desired properties of the fraction (for instance, the fact 
that a term is centered and/or some terms are mutually orthogonal, see Proposition 3 of 
[Pistone and Rogantin| | |2008[ )). 

In this a way the number of variables involved in the problem is reduced and the actual 
computation of all admissible patterns of a fraction can be simplified, as shown in the following 
examples. 

We consider orthogonal arrays of strength 2, i.e., with balanced simple and second order 
interactions. 

• Fractions of a 3^ design with 9 runs. 

The fractions have 9 runs. For fc = 0,1, 2, and a,b,c= 1, 2, we have: 

^ — 0 'k^{a.,h,c) ,k 9 , 

0 ’^( 0 , 0 , 0),0 = 9 : '>T'( 0 , 0 , 0),1 = '' 1 ( 0 , 0 , 0),2 = 0 , 

0 ^l(2,[3-a],[3-6]),[3-A:] = W(i,o,6),A: 

0 ^(a,0,0),fc 'O(^0,a,0),k ^(0,0,a), A: 3 , 'O^a,b,0),k ^(a,0,6),fc '^(0,a,b),k 3 

The relationships of Corollary[^ computed with CoCoA (see 
( [Abbott et al.[ )), are: 

^{i,i,i),o “ ''^^ 1 , 1 , 1 ),! ~ ~ 9n(i^i^i)_o + 9n(i,i^i)^i, 

~ 2 n(i^i^i)_o»^(i,i,i),i + + 2n(i_i^i)^o’^{i,i,i),2 ~ ''^(1,1,1),2 ~ 9 n(i^i^i)^i -|- 9 n(i_i^iy 2 

In order to easily handle these polynomials we compute the Grobner basis of corresponding 
ideal: 

^( 1 . 2 , 2),0 + ^( 1 , 2 , 2),1 + ’ 1 -{ 1 , 2 , 2),2 ~ 9, W(i_ 2 ,l ),0 + ’^( 1 . 2 , 1 ),! + ’T-(l, 2 ,l ),2 ~ 9, 

'T-( 1 , 1 , 2),0 + ^( 1 , 1 , 2),1 + ? 1 { 1 , 1 , 2),2 ~ 9, ’T-( 1 , 1 , 1),0 + ~ 9, 

^( 1 , 1 , 1),2 3” 9n'(i^i^i)^25 ^( 1 , 1 , 1),2 ^^^( 1 , 1 , 1),2 “b 2771(1^1^1)^2; 

’^(1,1,1),1* *^(1,1,1),2 + 1/2?J'(i,i,i), 2 “ 9/2n(i_i_i)^2 ■ 


Abbott, Bigatti, and Lagorio 










1.4. GENERALIZED WORD LENGTH PATTERN 


The admissible configurations for are: 

(3,3,3), (9,0,0), (0,9,0), (0,0,9). 

The counts for the other a follow from the constraints C1-C4 listed above. 

• Fractions of a 5^ design with 25 runs. In this case the computation is heavier than in 
the previous example. The constraints C1-C4 give 32 polynomials in 80 variables (the 
counts), one of which is shown below. 
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”( 1 , 4 , 3),3 “ 

”( 1 , 2 , 1),2 

”( 1 , 4 , 3 ) 

,4 

+ 

”( 1 , 2 , 1),4 

”( 1 , 4 , 3),4 “ 

”{ 1 , 3 , 2),3 

”( 1 , 4 , 3),4 

+ 

”( 1 , 3 , 2),4 

”( 1 , 4 , 3),4 

+ ”( 1 , 1 , 1),0 

”( 1 , 4 , 4),0 “ 

”( 1 , 1 , 1),2 

”( 1 , 4 , 4 ) 

,0 

+ 

”( 1 , 2 , 2),0 

”( 1 , 4 , 4),0 “ 

”{ 1 , 2 , 2),3 

”( 1 , 4 , 4),0 

+ 

”( 1 , 3 , 3),0 

”( 1 , 4 , 4),0 

“”( 1 , 3 , 3),4 

”( 1 , 4 , 4),0 + 

”( 1 , 1 , 1),1 

”( 1 , 4 , 4 ) 

,1 

- 

”( 1 , 1 , 1).3 

”( 1 , 4 , 4),1 + 

”( 1 , 2 , 2),1 

”( 1 , 4 , 4),1 

- 

”( 1 , 2 , 2),4 

”( 1 , 4 , 4),1 

“”( 1 , 3 , 3),0 

”( 1 , 4 , 4),1 + 

”■( 1 , 3 , 3 ),! 

”( 1 , 4 , 4 ) 

,1 

+ 

”( 1 , 1 , 1),2 

”( 1 , 4 , 4),2 “ 

”( 1 , 1 , 1),4 

”( 1 , 4 , 4),2 

“ 

”( 1 , 2 , 2),0 

”( 1 , 4 , 4),2 

+ ”( 1 , 2 , 2),2 

”( 1 , 4 , 4),2 “ 

”( 1 , 3 , 3),1 

”( 1 , 4 , 4 ) 

,2 

+ 

”( 1 , 3 , 3),2 

”( 1 , 4 , 4),2 “ 

”{ 1 , 1 , 1),0 

”( 1 , 4 , 4),3 

+ 

”( 1 , 1 , 1),3 

”( 1 , 4 , 4),3 

“”( 1 , 2 , 2),1 

”( 1 , 4 , 4),3 + 

”( 1 , 2 , 2),3 

”( 1 , 4 , 4 ) 

,3 

“ 

”( 1 , 3 , 3),2 

”( 1 , 4 , 4),3 + 

”( 1 , 3 , 3),3 

”( 1 , 4 , 4),3 

“ 

”( 1 , 1 , 1),1 

”( 1 , 4 , 4),4 

+ ”( 1 , 1 , 1),4 

”( 1 , 4 , 4),4 “ 

”( 1 , 2 , 2),2 

”( 1 , 4 , 4 ) 

,4 

+ 

”( 1 , 2 , 2),4 

”( 1 , 4 , 4),4 “ 

”( 1 , 3 , 3),3 

”( 1 , 4 , 4),4 

+ 

”( 1 , 3 , 3),4 

”( 1 , 4 , 4),4 


In this example it is unfeasible to explicitly compute all the admissible configurations, 
nevertheless one can check whether a given set of counts is admissible or not. 


1.4 Generalized Word Length Pattern 


As mentioned in the Introduction, the j-ih element of the GWLP of a fraction can be computed 
as the sum of the aberrations of the interaction terms A“ of order j. In this section, we make use 
of the formula in Equation 0 to express the aberrations as functions of the level counts, and 
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we study the mean aberration of a term over the level permutations. The analysis of the mean 
aberration is valid for generic asymmetric multilevel designs, while its use for the computation 
of the GWLP is limited to symmetric s’" designs, with s prime number. 

Definition 3. Given an interaction defined on a fraction T of the full factorial design D, 
its aberration, or degree of aliasing, Oq is given by the real number 

II., I|2 

II II 2 


where ||a ;||2 is square of the norm of the complex number x. 

The GWLP A{iF) = (^i(7^),..., of a fraction T is defined as 


aa:f)= Y. 


j = l,...,m 


ll“llo=J 

where ||a||o is the number of non-null elements of a, i.e., the order of interaction of 

The following proposition allows us to compute the aberration Uc without using complex 
computation. 

Proposition 4. Let be a simple or interaction term with values in fit. Its aberration Oa is 


1 


E- 


271 , 


k j Y, 




Proof. This formula follows from Dehnition [fusing Equation 
We observe that, if t is odd, we get 
1 


(1.4) 

□ 


G-i {«-i)/2 

E<« + 2 E cos 

I 7—0 k—1 


kj ^ ^ riQ,,7n.Q,jt—fc] 

2 7=0 


Moreover, for t = 2, 3, 4 we obtain 


n a„ = < 


Uq + Ui — 2noni 


if f = 2 


no -I- rii -b 712 — 720712 — nino — n2ni if f = 3 
Uq + nf + n 2 + — 2non2 — 2nins if f = 4 


where the suffix a is omitted to simplify the notation. 

It is well known that if a term is balanced on a fraction (i.e., Ua.k = ri/t for all 


k = 1,... ,t — V), then Oc = 0. Equation (1.41 gives an alternative proof of this fact. However, 


two disadvantages in the definition of Oa come out when facing with qualitative factors: (a) the 
aberration Oa is not independent on the permutation of the levels; (b) the fact that Uq = 0 does 
not guarantee that term with t levels, is balanced if t is not prime. 

To illustrate such problems, let us consider a term X°‘ with 6 levels and counts Ua. = 
(uo -b h, ui,uo, ui + h, uq, ui) for given uq, Ui, h € N. It is straightforward to check that = 0, 






1.4. GENERALIZED WORD LENGTH PATTERN 


Ui + h 



Figure 1: Three permutations of counts for a 6-level term. Left: Ua = 0; center: Uq, = 
{{uq — Ml) + h)'^jn?] right: jri^ 

even though is not balanced. Moreover, permuting the first two counts (ui ,mq + h,uo,u-i + 
h,uo,ui), Equation IH} yields n^aa = ((uq — ui) + h)^, while permuting the first and the 
third counts, (uo,Ui,uo + h,ui + h,uo,ui), one obtains = Sh^. To show the relevance of 
the position of the counts, such three configurations of counts are depicted on the unit circle in 
Figure |1.4| 

1.4.1 Mean aberration 

As a consequence of the previous discussion, we define here a permutation-invariant aberration 
of a term A“ as the mean of the aberration obtained through all the possible permutations of 
the counts {no,..., nt_i} of 

Definition 4. Given a fraction oi T> let A“ be a simple or interaction term with t levels 
and be {riafi, ■ ■ ■, the number of occurrences of its levels (0,1,..., t — 1) respectively. The 

mean aberration of is the mean of all the a^’s obtained by permuting (n^.o, • • ■ 


1 



(1.5) 


where 7r(na,o,.. •, na,t-i) denotes all the permutations of the counts of A“. 


The subscript a in na,i will be omitted in the following when unnecessary. 


Proposition 5. The mean aberration is: 



( 1 . 6 ) 


Proof. Take Equation (1.4| and first separate the cases k — 0 and k > 0 obtaining: 
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The first addendum is permutation-invariant and therefore it is enough to rewrite the second 
addendum. For fixed fc 7 ^ 0, we get 

t-i 

‘njnii-k] 


1 

ti 


n(no,...,nt_i) i = 0 

where in the summations there are t(t\) addenda. For fixed i and j each monomial niUj with 
i ^ j appears exactly t{t — 1 ) times, so that 


t-i 

7r(no,...,nt_i) i=0 


t-1 t-1 

E 


t -1 


nitij 


j7^i,j = l 


The first equality in Equation (1.51 is then proved, recalling that X]I=i ~ 

The second equality is straightforward. □ 

Remark 1. Notice that if the number of levels of X°‘ in I? is a prime number, the mean 
aberration of X°‘ on a fraction is zero if and only the corresponding aberration is zero: 

Oo, = 0 flc = 0 no = ni = ■ • ■ = rit-i ■ 

It is an easy consequence of the properties of the roots of the unity. On the other hand, if the 
number of levels is not a prime number, then aberration and mean aberration have different 
behaviors. The mean aberration is zero if and only if all the levels appear equally often: 


Ua = 0 no = ni = ■ ■ ■ = nt-i 

see the right expression of Equation |1.6| Whereas this fact is not true for aberration, as explained 
in examples of Figure [T^ where t = 6 and the level counts are not equal. In all the three cases 


Uc = ^^(9 (mo - wi) + 8/1 ) 


while the three aberration are 


Uq, = 0 , tta = ((mo — Ml) + h)^ fn^ , ac = ?>]%' jr? . 

Another nice property of mean aberration is that it represents the variances of the counts 
{no,..., ns_i} as shown in the following proposition. 

Proposition 6. Given a fraction T of T) let be a simple or interaction term with t levels 
and with level counts {no,..., nt-i}. The mean aberration Oa of X°' is the variance of 
{no,... ,nt_i}. 

Proof. From Proposition!^ we get 


1 1 


n^ t 




i=0 j=i 

t-1 


i=0 j=Q 








\2 

n, = 


n 

¥ 


^{t-1) 
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y^t — 1 ^2 ^2 

where a = —^ variance of {no,..., nt-i}. 


□ 


The variance is an index to describe the distribution of counts (no, ni,..., rit-i). Other 


indices can be used and, for instance, the authors in Fontana et al. 12014c I shown the connections 


between the aberration and Gini concentration index in the two-factor multi-level case. 

Remark 2 (Class of equivalence of aberrations). The value of Uq does not change for permuta¬ 
tions of all the counts in conjugate positions, e.g., for t = 5, permuting ni with 714 and 712 with 


773 simultaneously. It follows from Equation (1.41. Moreover the same holds for each of the t — 1 


circular permutations of levels, e.g., for t = 5, permuting (tiq, Tii, 712,773, 714) in (ni, 772,773,774, no) 
or (n 2 , ns, n 4 , no, ni) or (ns, n 4 , no, ni, ns) or (n 4 , no, ni, ns, ns). Then, there are {t— l)!/2 class 
of permutation with different Oq. 

We explicitly show that the permutation of only two counts in conjugate position is not 
enough to preserve the same aberration. Consider two terms Af and X 2 with t = 5 levels 
and counts (no, ni, ns, ns, n 4 ) and (no, n 4 , ns, ns, ni) respectively, where only ni and n 4 are 
permuted. Let and the two corresponding aberrations. The difference between al, and 


27r 477 

2 I cos -cos — ) (ni — n4)(n2 — ns) . 


1.4.2 GWLP 

From now on we consider symmetric fractional factorial designs s’" with s prime, where each 
A“ takes values in 

Restricting to symmetric designs with a prime number of levels, GWLP of a fraction A, 
Ajr = (Ai(A),..., Am{J-)), can be computed using only the mean aberrations. 


Proposition 7. Let T be a fraction of a s’" design, s prime number. Then j-th element of the 
GWLP is: 


ldlo=J 


Proof. Let A“ be a term of order j with counts (no,... ,ns-i). Consider the powers of X°‘ 
with a multi-exponents in L, namely X°‘, ..., Observe that, as s is prime, the 

counts of such powers are all permutations of (no,..., n^-i). In more details, these counts have 
the following properties. 


• no is always in the count associated to ujq in all the powers h £ {1 ,..., s — 1}. 


• Let ni, 7 7^ 0, be the count associated to uji of X°‘. m is the count associated to uJ[hi\ of 
and to ui[ki] of for h ^ k £ {1,..., s — 1}. Since s is prime, this means that, 

for each i, ni becomes the count associated to a different root of the unity exactly once 
for all the powers. 
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Therefore, the sum of the aberrations of the s terms above rewrites as: 

cia + a[2a] + • • • + a[(s_l)a] = ^ ^ \ j , 

where the permutations tt range over the s permutations defined above. As in the proof of 
Proposition!^ the first addendum is permutation-invariant and therefore it suffices to rewrite 
the second addendum. Now notice that in the sum 


s-l 

7r(no , • • • ,^s — 1) ^—0 

for all fe 7 ^^ 1 each monomial riirij with i ^ j appears exactly once, so that 


O-a + a[2Q:] + • • • + U[(s —l)a] — ~ l)<Ia • 


This completes the proof. □ 

Remark 3. We emphasize again that the advantage of using mean aberrations in place of the 
classical aberration formula is the computational ease. In fact, it is enough to compute the s’" 
vectors of counts, one of each term x°‘, a £ L, and then use the terse formula in Equation ( |1.6[ ). 
Moreover, according to the previous discussion, the s’" mean aberrations are equal s by s and 
this reduces further the computational cost. 

Remark 4. In the non prime case the GWLP is not the sum of the mean aberrations, as shown 
in the following counter-example. Consider a trivial example, namely a design with only one 
factor X with four levels and, on a fraction with 6 runs ni = (1, 2,1, 2). Then has 2 levels 
with n 2 = (2,4), and X^ has 4 levels with na = (1, 2,1, 2). We have: 

ai = as = 0, 02 = 1/9, 01=03 = 1/27, 02 = 1/9 . 

Then Ai — ~ th® sum of the mean aberrations is 5/27. 


1.5 Case studies 


In this section we present some results concerning the use of mean aberrations to distinguish 
among different fractions. We mainly focus on the most interesting situation where fractions 
have the same GWLP. We use some non-isomorphic orthogonal arrays extracted from the com¬ 


plete series of non-isomorphic orthogonal arrays, Schoen et al. (20101. We consider examples of 


orthogonal arrays of strength 2 , so that Ai[X) = A 2 {J-) = 0 in all cases, and the first interesting 
aberrations correspond to the terms of order 3. 






1.5. CASE STUDIES 


1.5.1 2-level orthogonal arrays, OA{16]2]2^^) 

There are six non-isomorphic orthogonal arrays of strength 2 with ten 2-level factors and 16 
runs. We denote these arrays by i = 1,..., 6. They have the same GWLP: 

(0,0,8,18,16,8,8,5,0,0). 

For each orthogonal array we compute the 10 mean aberrations of order 3 and the 210 aberra¬ 
tions of order 4. The mean aberrations corresponding to the interaction terms of order 3 can 
be classified as in Table [T] 

Table 1: Aberrations of the terms of order 3 for the OA(16; 2; example. 


^OL 

.Ai 

:f2 

.A 3 


A-s 

As 

0 

112 

100 

100 

88 

88 

88 

0.25 

0 

16 

16 

32 

32 

32 

1 

8 

4 

4 

0 

0 

0 


Table shows that the six orthogonal arrays can be clustered into three groups: J-i, 
and {7^4, The mean aberrations corresponding to the interaction terms of 

order 4 can be classified as in Table Here the six fractions are now six different distributions. 
Moreover, one can see that contains a regular fraction with 8 defining words of order 3 and 
18 words of order 4. Finally, we remark that the two values of the GWLP = 8 and 

A 4 ,{Fi) = 18 may be recovered from the mean aberrations via a simple weighted sum. 

Table 2: Aberrations of the terms of order 4 for the OA(16; 2; 2^°) example. 



Ai 

A2 

A3 

A4 

As 

Ae 

0 

192 

168 

180 

192 

168 

180 

0.25 

0 

32 

16 

0 

32 

16 

1 

18 

10 

14 

18 

10 

14 


1.5.2 3-level orthogonal arrays, 074(18; 2; 3^) 

There are three non-isomorphic orthogonal arrays of strength 2 with seven 3-level factors and 
18 runs. We denote these arrays by Fi, F 2 and F^. They have the same GWLP: 

(0,0,22,34.5,27,31,6). 

For each orthogonal array we compute all the 278 values of the mean aberrations to the inter¬ 
action terms of order 3. They can be classified as in Table In this case the distribution of the 
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mean aberrations corresponding to the interaction terms of order 3 is sufficient to distinguish 
the arrays. 

Table 3: Aberrations of the terms of order 3 for the OA(18; 2; 3^) example. 


1 

J-i 

^2 

^3 

0 

134 

198 

102 

0.083 

96 

0 

144 

0.25 

48 

80 

32 


1.5.3 5-level orthogonal arrays, 0^4(25; 2; 5^) 

There are two non-isomorphic orthogonal arrays of strength 2 with three 5-level factors and 25 
runs. We denote these arrays by T\ and J- 2 - They have the same GWLP: 

(0,0,4). 

For each orthogonal array we compute the 64 values of the mean aberrations corresponding 
to the interaction terms of order 3. They can be classified as in Table Also in this case 
the distribution of the mean aberrations corresponding to the interaction terms of order 3 is 
sufficient to distinguish the arrays. 

Table 4: Aberrations of the terms of order 3 for the OA(25; 2; 5^) example. 


Q-ct 



0 

12 

60 

0.04 

16 

0 

0.06 

32 

0 

0.36 

4 

0 

1 

0 

4 


1.6 Future directions 

The theory presented in this paper and the examples discussed in the previous section show 
the relevance of the notion of mean aberration to discriminate between designs with the same 
GWLP, but several problems are still open. We mention here only a couple of directions for 
future work. First, it should be interesting to generalize the notion of mean aberration in 
order to apply it in the general case, namely for factors with number of levels not prime, and 
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for asymmetric designs. Second, we deem important to explore the connections between the 
geometric structure of the design points of the fraction and its mean aberration. In fact, some 
results in this direction have already been achieved, and they are presented in [Fontana et alT] 
( |2014a[ ), [Fontana et al.| | [2014b[ ), and [Fontana et ^ ( |2015| . In those papers, saturated fractions 
and D-optimal saturated fractions are described in terms of the circuit basis, a combinatorial 
object computed from the model matrix. Although aberration and GWLP are defined in a 
model-free framework, nevertheless we think that the study of the geometry of the fractions will 
yield new interesting results on aberration, and new insights in design selection. 
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